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The charm fragmentation function has been measured in photoproduction of jets con- 
taining D*^ mesons. The measured function has been used to extract free parameters 
of different fragmentation models. Measurements of excited charm, _Di(2420)'' and 
-D5(2460)", and charm-strange, 1)^(2536)=^, mesons and a search for the radially ex- 
cited charm meson, _D*'(2640)*, were also performed. The results are compared with 
those measured previously and with theoretical expectations [T]. 

1 Introduction 

Charm hadrons were produced copiously in ep collisions with a centrc-of-mass energy of 
318 GeV at HERA providing a means to study charm hadronisation. During first phase of 
the HERA operation (1992-2000), the ZEUS collaboration accumulated data sample corre- 
sponding to an integrated luminosity of ~ 120 pb~^. Measurements of the and 
A+ productior0 with the HERA I data were used to determine charm fragmentation ratios 
and fractions of c quarks hadronising as a particular charm hadron, /(c — > D, Ac), in earlier 
ZEUS studies |H [3] . Recent ZEUS measurements of the charm fragmentation function [i] 
and the production of excited charm and charm-strange mesons [5] are summarised in this 
note. 

2 Measurement of charm fragmentation function 

The measurement of the charm fragmentation function in the transition from a charm 
quark to a £)*"*" meson was performed in photoproduction regime with the virtuality of 
the exchanged photon < 1 GeV^ and the photon-proton centre-of-mass energy in the 
range 130 < W < 280 GeV. The D*~^ mesons were reconstructed from the decay chain 
£)*+ I?°7r+ — > {K^TT^)^^ using the mass difference technique. The D*~^ meson was 
included in the jet-finding procedure and was thereby uniquely associated with one jet only. 
The fragmentation variable, z, was defined as 

(£;+P||)^'V(i?+P||r' ^ {E+p\\)^'^/2Ei'^\ 

where p\\ is the longitudinal momentum of the D*~^ meson or of the jet relative to the axis of 
the jet of energy E^'^*^. The equivalence of {E +p\\y^^ and 2 E^°^ arises because the jets were 
reconstructed as massless objects. The measurement of the normalised differential cross 
section, l/a{da/dz), was performed in the kinematic range pt{D*^) > 2 GeV, \r]{D*^)\ < 
1.5, i;^* > 9 GeV and {rf'^^l < 2.4. The above requirements on pt{D*+) and i;^* allowed 
the fragmentation function measurement in the range 0.16 < z < 1. 

* Participation in DIS 2009 was supported by the Deutsches Elektronen-Synchrotron DESY 
^Hereafter, charge conjugation is implied. 
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Figure 1: Charm fragmentation function in transition to D*"*" for the ZEUS data compared 



to (left) measurements in 



annihilations and (right) predictions of FMNRxCj 
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The measured fragmentation function is compared in Fig. [TJ^left) with previous mea- 
surements from Belle [6], CLEO [7] and ALEPH [8j. For shape comparison, the data sets 
were normalised to l/(bin width) for z > 0.3. The Belle and CLEO data are measured at 
a centre-of-mass energy of ^ 10.5 GeV, whereas the ALEPH data were taken at 91.2 GeV. 
The corresponding scale of the ZEUS data is given by twice the average transverse energy 
of the jet, 23.6 GeV. The ZEUS data in Fig. [ijleft) are shifted somewhat to lower values 
of z compared to the CLEO and BELLE data with the ALEPH data even lower that is 
consistent with expectations from scaling violations in QCD. 

The ZEUS data were compared with fragmentation models implemented in the leading- 
logarithmic Monte Carlo (MC) program PYTHIA [9j. The LUND string fragmentation 
model [To] modified for heavy quarks [TT] gives a reasonable description of the data. The 
PYTHIA predictions obtained using the Peterson fragmentation function [T^] was fit to the 
data via a x^-minimisation procedure to determine the best value of the parameter e. The 
result of the fit is e = 0.062 ± 0.0071qoq4. The result is in reasonable agreement with the 
default value used in PYTHIA (0.05), with the value measured by the HI collaboration 
in deep inelastic scattering (0.061!I!q qoq) [13], and with the value 0.053 obtained in the 
leading- logarithmic fit [H] to the ARGUS data [15] . 

The data were also compared with the next-to-leading-ordcr (NLO) QCD predictions 
from Frixionc ct al. (FMNR) [16]. The predictions with the parton-level jets were translated 
to the predictions with the hadron-level jets using the hadronisation correction factors, 
^had"""' obtained with the PYTHIA MC. The result of varying e in the Peterson hmction 
for the predictions of FMNRxC^g^""" is shown in the Fig. [TJright). The fit result, e = 
0.079 ± 0.008j:!^:^J^, exceeds the value 0.035 obtained from the NLO fit [M] to the ARGUS 
data [15]. The fit of the FMNRxCjg^"'" predictions with the fragmentation function from 
Kartvelishvili et al. [17] yielded the value of the free parameter a = 2.67 ± O-lSj^Q ^s- 
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3 Production of excited charm and charm-strange mesons 



The first study of excited charm and charm-strange mesons produced in ep coUisions was 
restricted to decays, for which significant signals were identified: 







"71" , 








+ 







The measurement was performed in the fuU kinematic range of Q^. 
were identified using the decay channel D*+ — > D'^tt^ following by either A'~7r+ or 

D° — > /^"Tr+Tr+Tr" decay. The Z?+ mesons were reconstructed in the decay /v "7r+7r+. 

The D*" mesons were tagged in the decay to a Z?" and undetected neutrals following by the 
D° K-TI+ decay. 

To extract the and -Dj" yields and properties, a minimal fit was performed using 
simultaneously the M(Z)+7r") distribution and the M(D*+7r") distributions in four helicity 
intervals. The helicity angle [a) is defined as the angle between the momenta of the addi- 
tional pion and the pion from the D*'^ decay in the D*'^ rest frame. The helicity angular 
distribution can be parametrised as 

oc 1 + /i cos a, (1) 



dcos 

where h is the helicity parameter. Heavy Quark Effective Theory [TH] (HQET) predicts 
/i = 3 (ft. = 0) for the l"*" state from the j = 3/2 (j = 1/2) doublet, and /i = — 1 for the 
2+ state from the j = 3/2 doublet. Only D-wave decays are allowed for the members of 
the j = 3/2 doublet; therefore they are supposed to be narrow. On the other hand, the 
members of the j = 1/2 doublet decay through S'-wave only and therefore are expected to 
be broader |19| . Due to the finite charm quark mass a separation of the two doublets is only 
an approximation and amplitudes of two observable states with = 1+ can be mixtures 
of D- and S'-wave amplitudes. 

The measured masses of the 1)° and f-re in reasonable agreement with the world 
average values [10] • The measured width is TiD^I) = 53.2 ± 7.2^4:9 MeV which is above 
the world average value 20.4± 1.7 MeV [5D]. The observed difference can be a consequence of 
differing production environments. The width can have a sizeable contribution from the 
broad S'-wave decay even if the S-wave admixture is small |21| . A larger S'-wave admixture at 
ZEUS with respect to that in measurements with restricted phase space, which can suppress 
production of the broad state, could explain why the measured width is larger than the 
world average value. 

The measured Dj* helicity parameter is h{Di) ~ S-Qli'yii'o which is inconsistent with 
the prediction for a pure S'-wave decay of the 1+ state, h — Q. It is consistent with the 
prediction for a pure D-wave decay, h = 3. In the general case of D- and S'-wave mixing, 
the helicity angular distribution form of the 1+ state is: 

dN 



dcosa 



ccr+{l- r)(l + 3cos^Q!)/2 -t- \/2r{l - r) cos 0(1 - Scos^a), (2) 



where r ~ '^s/ij's + ^d), Ts/d is the S- / D-wave partial width and (j) is the relative phase 
between the two amplitudes. Using Eqs. ([T]) and cos0 can be expressed in terms of r 
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Figure 2: Cosine of the relative phase of S- and D-wave amphtudes versus r = r5/(r5 + ro) 
in (left) the Di(2420)" ^ D*+Tr- decay and (right) the D^i(2536)+ ^ D*+K'^ decay. 



and the measured value of the hclicity parameter, h: 



{3-h)/{3 + h)~r 

coso— = . (S) 

2v/2r(l-r) 

Figure [2{left) compares with previous measurements the range restricted by the measured 
h{Di) value and its uncertainties in a plot of cos (j) versus r. The ZEUS range has a marginal 
overlap with that restricted by the CLEO measurement of h{D'() = 2.74lo 9° [22]. BELLE 
performed a three-angle analysis and measured both the cos (f> and r values [22 . The BELLE 
measurement, which suggested a very small admixture of .S-wave to the Di (2420)" D*^n~ 
decay and almost zero phase between two amplitudes, is outside the ZEUS range; the dif- 
ference between the two measurements, evaluated with Eq. ([3]), is 2 standard deviations. 

A signal from the D*'^K^ decay was observed in the Af(£'"A'+) distribution 

with an average negative shift of 142.4 ± 0.2 McV with respect to the nominal mass. 
To extract the D^i yields and properties, an unbinncd likelihood fit was performed using 
simultaneously values of M{D^K^), M {D*^K^), and cos(q!) for Z)*+A'^ combinations, with 
the helicity angle defined as the angle between the momenta of A'^ and the pion from the 
decay in the D*'^ rest frame. 

The measured Df^ mass is in good agreement with the world average values [20]. The 
measured Df^ helicity parameter is h{D'^^) = ^^■7'i-oiijlLo''o5- "^^^ measured h value is 
inconsistent with the prediction for a pure D-wave decay of the 1+ state, /i = 3, and is 
barely consistent with the prediction for a pure S-wave decay, h — 0. Figure Upright) shows 
a range, restricted by the measured h{D'^^) value and its uncertainties, in a plot of cos0 
versus r = rs/(rs -|- F^j) (Eq. |3|). The measurement suggests a significant contribution 
of both D- and S-wave amplitudes to the £>si(2536)+ D*+ Kl decay. The ZEUS range 
agrees with that restricted by the CLEO measurement of h(Df-^^) = — 0.23^0 32 [21] ^-nd with 
the BELLE three-angle measurement of both cos (f> and r values [25] . 

The measured yields were converted to the fragmentation fractions f{c^ Dl) = 3.5 ± 
0.4t°;^%, /(c ^ Df) = 3.8 ± G.7l°J% and /(c ^ D+) = 1.11 ± G.16l°°o%- The 

annihilations. The measured ratios of 



DIS2009 



the dominant ^^.d D^j^ branching fractions are 

— flJ^ = 2.8 ± 0.8t°-^ , =2.3 ±0.6 ±0.3 

in agreement with the world average values [20| . 

No radially excited meson, reported by DELPHI [26], was observed. An upper 

limit, stronger than that obtained by OPAL [57], was set on the product of the fraction of 
c quarks hadronising as a £)*'+ meson and the branching fraction of the L)*+7r+7r^ 
decay in the range of the mass from 2.59 to 2.69 GeV: 

/(C ^ D*'+) ■ BD-'+^D'^rr^rr- < 0.4% (95% C.L.). 
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